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The construction of the operators and correlators required to determine the excited baryon spectrum is pre- 
sented, with the aim of exploring the spatial and spin structure of the states while minimizing the number of 
propagator inversions. The method used to construct operators that transform irreducibly under the symmetries 
of the lattice is detailed, and the properties of example operators are studied using domain-wall fermion valence 
propagators computed on MILC asqtad dynamical lattices. 



1. INTRODUCTION 

The determination of the excited baryon spec- 
trum is an important step on the road to under- 
standing the dynamics of QCD, and it is, for ex- 
ample, a vital element of the Jefferson Laboratory 
experimental program. Lattice gauge theory has 
a pivotal role not only in providing ab initio cal- 
culations of the masses of the lowest lying states, 
but also in establishing the quark and gluon con- 
tent of the states. There has thus been increas- 
ing activity amongst the lattice community that 
has demonstrated that the masses of the lowest- 
lying baryon states of both parities for spin- 1/2 
and spin-3/2 are accessible to lattice calculation, 
and there have been calculations of the mass of 
the first positive-parity, spin-1/2 excitation of the 
nucleon, the so-called Roper resonance. 

Most of these calculations have been made in 
the quenched approximation to QCD, and an im- 
portant component of future lattice studies will 
be a careful analysis of the systematic uncertain- 
ties on these results. However, the confrontation 
of lattice calculations with experiment will re- 
quire a more complete understanding of the spec- 
trum, obtaining both the masses of the higher 



spins and of the radial excitations. Furthermore, 
the discovery of the quark and gluon structure of 
the excited states will require a more extensive 
basis of interpolating operators than used so far, 
allowing both for excited glue and for multi-quark 
operators. The aim of this talk is to describe the 
design of such operators, and in particular to con- 
struct operators that transform irreducibly under 
the symmetries of the lattice. 

The remainder of this talk is laid out as fol- 
lows. In the next section, the methodology of 
computing the spectrum is described, and the 
symmetries used to classify baryon states are de- 
tailed. Gauge-invariant "elemental" operators 
that have the correct flavor structure and that 
explore the spatial structure of the baryons are 
then introduced. The reduction of these opera- 
tors to ones that transform irreducibly under the 
cubic symmetry of the lattice is then performed. 
Finally, the reduction is illustrated for the case 
of point-like elemental operators obtained using 
domain-wall fermion (DWF) valence propagators 
computed on a dynamical asqtad-fermion back- 
ground. 
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2. THE SPECTRUM FROM LATTICE 
CALCULATIONS 

The computation of the spectrum of states in 
lattice QCD is in principle straightforward. 

1. Choose an interpolating operator O that 
has a good overlap with P, the state of in- 
terest, 

(0 I O | P) £ 0, 

and ideally a small overlap with other states 
having the same quantum numbers. 

2. Form the time-sliced correlation function 



3. Insert a complete set of states between O 
and O*. The time-sliced sum puts the inter- 
mediate states at definite momentum, and 
we find 
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where the sum over P includes the contri- 
butions from two-particle and higher states. 

4. Continue to Euclidean space t — > it, yield- 
ing 



Cit) = £ 
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so that the correlator falls off exponentially 
with the mass of the lightest state at large 
times. 

In order to obtain the mass of the higher reso- 
nances in Eqn.^ the most appealing approach is 
to apply variational methods [02] , beginning with 
the computation of a matrix of correlators 



(2) 



where the {Oi} form a basis of operators of defi- 
nite quantum numbers. In the usual application 
of the variational technique, we find a basis of 
eigenvectors V(to) that diagonalizes the transfer 
matrix C(to)~ 1 C(t) for t = to+1 for some to close 
to the source, and determine the spectrum from 
the eigenvalues of 



v(t + ir'cito^c^vito + 1) 



(3) 



at subsequent times. The method relies on con- 
structing a basis of operators {Oi} that provides 
a good description of the states of interest, and 
the remainder of this talk will focus on the con- 
struction of that basis. 

We classify the states of the continuum by their 
flavor structure, F, and by their parity P, total 
spin J and "helicity" J z . Whilst the flavor struc- 
ture can be faithfully preserved in a lattice cal- 
culation, the parity, spin and helicity labels arise 
from looking at the symmetry properties of states 
under rotations and reflections. In a lattice calcu- 
lation, we have replaced continuum space time by 
a hypercubic lattice, and thus must classify states 
according to the irreducible representations A of 
the cubic group (for particles at rest); the row A 
within the representation is then the analogue of 
helicity. 

In the following section, we will construct a set 
of elemental baryon operators {Bf(x)} having 
the correct flavor properties F. We will then use 
group theory to generate a set of baryon operators 
{E>[ AX (x)} that transform irreducibly under the 
rotation and reflection symmetries of the lattice, 
and from which we will construct our correlation 
matrix 



C^(t) 



A ^(x)Bf F (0)), 



(4) 



and hence extract the spectrum of states corre- 
sponding to A, A. 

3. ELEMENTAL BARYON OPERA- 
TORS 

The color and flavor structure of the operators 
is dictated by the requirements of gauge invari- 
ance and of isospin respectively; it is natural in a 
lattice calculation to assume exact isospin sym- 
metry, and to classify states and operators by 
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their isospin and strangeness, rather than accord- 
ing to SU (3) flavor. 

The use of "smearing" has proven essential to 
reduce the coupling of operators to the higher 
excited states, and thus to ensure that the cor- 
relators Eqn. [5] are dominated by only a few 
states even close to the source. In the example 
given here, we will adopt a gauge-covariant spa- 
tial Laplacian to smear the fields: 

4>(x) = (1 + a 2 A/4N) N ip(x) (5) 

with the three-dimensional Laplacian defined by 

A^(x) = J2 (u k {x)^(x+k) - V>(a:)) , (6) 

k=±l,±2,±3 

where U denotes a link variable that can 
be smeared according to, say, the APE 
prescription 2|. Both TV and a are tunable param- 
eters, and in the limit N — > oo, Eqn. reduces 
to Gaussian smearing of width a. Note that the 
square of the smeared field vanishes, in the man- 
ner of a simple Grassmann field. 

An important consideration is to construct op- 
erators that enable us to build up the radial and 
orbital structure of the states. This is accom- 
plished by additional powers of the Laplacian, 
and by spatial derivatives, respectively, so that 
the quark building blocks become: 

X n Aaa (x) = (A>(x)) , (7) 

V / Aaa 

CIM = (bf^ix)) , (8) 

\ J / Aaa 

where A is a flavor index, a is a color index, a a 
Dirac spin index, and n is an integer. The p-link 
gauge-covariant forward displacement operator is 
defined by 

bfo{x) = 

U s (x) . . . U s {x + {p-\)j)0{x+ri), (9) 

where j — 1, 2, 3, and p is an integer specifying 
the length of the displacement. 

From these quark building blocks, we construct 
our elemental three-quark operators: 

VABC £ abc XAaa XBb/3 XCcyi 

(10) 

0ABC £ abc X^Aaa Xb\@ £cc7 ' 
VABC t abc XAaa ^Bb/3 ?Cc 7 • 



The flavor structure is specified by (fi^BC whilst 
the Levi-Civita symbol ensures that these are 
color singlets. The radial and orbital structures 
are explored by varying ni,ri2,ns, and through 
the choices of pi,P2>jik respectively. A large 
number of baryons can then be studied using 
a somewhat small number of quark propagator 
sources. 

4. SYMMETRIES OF THE LATTICE 

Proper rotations restricted to an isotropic cu- 
bic lattice form the cubic group O. It has 
24 elements, and five conjugacy classes and 
thus five single- valued irreducible representations: 
Ai,A2,E, T\,T2, of dimensions 1, 1, 2, 3 and 3 re- 
spectively. If we allow for spatial inversions, cor- 
responding to P = ±1, we obtain the group Oh, 
and the irreducible representations acquire a fur- 
ther label, g or u, corresponding to positive and 
negative parity states, respectively. The con- 
struction of operators transforming irreducibly 
according to the single-valued representations is 
well known, and has been crucial to identifying 
the states in lattice calculations of the glueball 
spectrum. 

There are four two-dimensional spinorial repre- 
sentations of Oh~. Gig, G±u, G2 g and G2 U , and two 
four-dimensional representatons, H g and H„ |4I5| . 
The irreducible representations J of the contin- 
uum group SU (2) are reducible under the cubic 
group O; the number of times that each of 
these reducible representations occurs in the irre- 
ducible representation r of O is shown in Tabled 
the extension to Oh is straightforward. 

Whilst full rotational symmetry is restored in 
the continuum limit, the irreducible representa- 
tions of Oh contain states of many representations 
of SU(2), and in general the degrees of freedom, 
corresponding to different hclicitics, lie in differ- 
ent irreducible representations. Thus, for exam- 
ple, a state of spin 5/2 has four degrees of freedom 
in the four-dimensional representation H, and 
two degrees of freedom in the two-dimensional 
representation G2', the identification of the spin 
of a state is accomplished by looking for the ap- 
proach to a common mass across the different ir- 
reducible representations in the continuum limit. 
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Table 1 

The number of times that the irreducible rep- 
resentation r of O occurs in the reduction of the 



irreducible representation J of SU(2). 


J 




n G 2 




1/2 


1 








3/2 








1 


5/2 





1 


1 


7/2 


1 


1 


1 


9/2 


1 





2 



Central to the task of finding the B^ XF (x) is 
the projection formula 



B ^f {S) = ^Y:d^\R)U r B({x)uI 
y° h Reo h 



(13) 



where A refers to an Oh irrep, A is the irrep row, 
go h is the number of elements in Oh, d\ is the 
dimension of the A irrep, Dml(R) is a A repre- 
sentation matrix corresponding to group element 
R, and Ur is the quantum operator which im- 
plements the symmetry operations; the temporal 
argument is suppressed. 

Application of this formula requires explicit 
representation matrices for every group element. 
Representation matrices for all allowed proper ro- 
tations can be generated from the matrices for 
C iy and Oi z , the rotations by n/2 about the y- 
and z-axes, respectively. The explicit matrices we 
use for these generators are given by 

1 



£> (Gl) (C 4 .) 



V2 
1 



2V2 



1 -1 
1 1 



l-i 
1+i 



D^(C 4y ), 
= -D^(C iz ), 



1 

3 



1 -a/3 V3 -1 

V3 -1-1 V3 

y/3 1 -1 -a/3 

1 V3 V3 1 

-l-i 

l-i 

1+i 

-1 + i 



5. EXAMPLE: LOCAL NUCLEON IN- 
TERPOLATING OPERATOR 

We will illustrate the procedure above by ex- 
amining the reduction of the local proton inter- 
polating operator, with / = 1/2, 13 = 1/2. The 
elemental operators that are gauge invariant and 
have the correct isospin are 



/p. abc 1 a ib c ia b c \ 



(14) 



Examination of Eqn. ^] reveals the constraints 
$ Q /3 7 + $ 7 /3q = and $ Q/ 3 7 + <&f 3la + $ 7Q(3 = 
so that there are only 20 independent operators. 
The projection operation Eqn. ^| is applied to 
each of the twenty linearly independent opera- 
tors in Eqn. 1141 and linearly independent opera- 
tors identified that transform irreducibly under 
Oh- The procedure depends on the Dirac ba- 
sis adopted, and therefore the reduction must 
be performed using the basis employed to com- 
pute the quark propagators. The reduction in 
the DcGrand- Rossi basis is shown in Table |2 the 
corresponding reduction in the Dirac basis is pro- 
vided in ref. [H]. 

There are three embeddings of G\ and a sin- 
gle embedding of H, of both parities. This is 
in accord with expectations from an analysis of 
the continuum operators, where it is possible to 
construct three linearly independent nucleon in- 
terpolating currents^, and a single spin-3/2 in- 
terpolating current. 

The computation of the correlator matrix 
Eqn.0]is straightforward. We compute the gauge- 
invariant "generalised" baryon correlator 

e abc e abc Q aa_( x . Q) Q t>b_^. Q)Q™( X] 0) (15) 

x 

where Q is the quark propagator, whence 
^{$ a /3y{x,t)$ & ^) = B a/3jA p v (16) 

X 

where B is obtained from B by symmetrizing over 
a and (3 to account for the two possible Wick 
contractions, and antisymmetrizing over a, 7 and 
over a, 7 to take account of the two terms in 
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Table 2 

Combinations of the operators * Ql a 7 in Ean. ll3| which transform irreducibly under Oh for the DeGrand- 
Rossi representation of the 7-matrices, employed by LHPC and MILC. 



Irrep 


Row 




Operators 




Gig 


1 




*112 + *334 




Gig 


2 




— *221 — *443 




Gig 


1 




*123 ~ *213 + *314 




Gig 


2 




*124 — *214 + *324 




Gig 


1 


2*114 + 


2*332 - *123 - *213 + 2*i 3 4 " 


" *314 


Gig 


2 


-2*223 " 


- 2*441 + *124 - *214 - 2*234 


+ *324 


Giu 


1 




*112 — *334 




Giu 


2 




-$221 + *443 




Giu 


1 




*123 — *213 — *314 




Giu 


2 




*124 — *214 — *324 




Giu 


1 


2*114 - 


2*332 — *123 — *213 _ 2*134 - 


f *314 


Giu 


2 


-2*223 - 


- 2*441 + *124 + *214 + 2*234 


_ $324 


H g 


1 




V3(*113 + *33l) 




Hy 


2 


*114 + 


*332 + *123 + *213 ~ 2*i 3 4 + 


*314 


H g 


3 


*223 — 


*441 + *124 + *214 - 2*234 + 


*324 


H g 


4 




\/3(*224 + *442) 




H u 


1 




\/3(*113 - *33l) 




H u 


2 


$114 - 


*332 + *123 + *213 + 2*i 3 4 ~ 


*314 


H u 


3 


*223 — 


*441 + *124 + *214 + 2*234 ~ 


*324 


H u 


4 




V3(*224 ~ *442) 





Eqn. Finally, we form the correlators 

^•(t) = ^rf^(*^ 7 (£,t)$^(o)), (17) 

X 

where the Ti , Tj are given in Table [21 While the 
computation of the generalized baryon correlators 
is computationally quite demanding, the subse- 
quent calculation of the 20 x 20 matrix of corre- 
lators of Ean. ll7| can be performed in a matter of 
minutes on a workstation, even for several hun- 
dred configurations. 

It is instructive to list some of the properties of 
the correlation matrix Eqn. 1171 

• Cross-correlations between operators 
within different irreducible representations 
vanish. 

• Cross-correlations between the operators of 
different rows within the same representa- 
tion vanish. 



• The different rows within a representation 
can be chosen to have the same normali- 
sation, so that correlators thus constructed 
have the same expectation values. 

• Cross-correlations between the same row of 
different embeddings of the same irreducible 
representation do not vanish; they can be 
constructed to vanish only on a particular 
time slice. 

Thus the strategy for extracting the mass spec- 
trum from the data is now clear. Firstly, choose 
a particular irreducible representation, possibly 
averaging over the rows within that representa- 
tion. Then form the matrix of correlators between 
all the averaged operators within that represen- 
tation, and extract the spectrum using, say, the 
method of Eqn. EH 

The talk concludes with an examination of the 
correlators on an ensemble of lattice data, and in 
particular a verification of the properties above. 
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To perform this study, we employ an ensemble of 
around 300 configurations computed on a 20 3 x 64 
lattice with Nf =2+1 flavors of asqtad staggered 
quarks having am u /am s = 0.01/0.050. The lat- 
tice spacing a = 0.124 fm is obtained from the 
IP-IS splitting in the Upsilon spectrum|U]. 

We employ domain-wall fermions at a single 
mass m — 0.01 for the valence quarks, which are 
computed on "chopped" 20 3 x 32 lattices, with 
N$ = 16; to reduce the residual mass, the configu- 
rations are first blocked using a non-perturbative 
HYP blocking scheme ^H]- Smeared sources are 
used, using a = 4.35 and N = 30 in Eqn. but 
local sinks, and thus the correlators are not pos- 
itive definite; full details of the calculation will 
appear in a subsequent paper jllj . 

The effective masses obtained from the diago- 
nal correlators using the three embeddings of G\ g 
and Giu are shown in figure^ The plot reveals a 
clear signal for both positive- and negative parity 
states from each of the three embeddings. Note 
that the three usual baryon interpolating oper- 
ators, for example in ref. correspond to lin- 
ear combinations of the operators used here; the 
identification of a "good" baryon operator corre- 
sponding to the nucleon, and a baryon operator 
having a good overlap onto the Roper resonance, 
should appear through the variational analysis. 
The effective masses for H g and H u are shown as 
the lower plot in figure ^ there is a clear sep- 
aration from the mass of the ground-state nu- 
cleon, with the negative-parity state being lower 
in mass, as observed experimentally. 

We illustrate the orthogonality properties in 
two ways in figure El firstly, by confirming that 
the cross-correlators between different irreducible 
representations are indeed zero, and secondly by 
verifying that those between different rows of the 
same irreducible representation vanish. The or- 
thogonality properties are clearly well satisfied in 
these representative examples, indicating that the 
numerical implementation is under control. 

6. CONCLUSIONS 

We have detailed the design of operators 
needed to determine the excited baryon spec- 
trum. The method used to construct operators 



transforming irrcducibly under the symmetries 
of the lattice has been outlined, and the effi- 
cacy of the method illustrated. The calculation 
of "generalised" baryon correlators followed by 
group-theory projections is computationally ef- 
ficient, and allows straightforward identification 
of the quantum numbers of the correlators. Fur- 
thermore, the methodology is readily extendible 
to more complicated operators, such as those for 
pentaquarks and those having excited glue. 

This work was supported by the U.S. Na- 
tional Science Foundation under Awards PHY- 
0099450 and PHY-0300065, and by the U.S. De- 
partment of Energy under contracts DE-AC05- 
84ER40150 and DE-FG02-93ER-40762. We are 
indebted to members of the MILC collaboration 
both for making available their configurations and 
for many valuable insights, and to Anna Hasen- 
fratz and John Negele for many helpful discus- 
sions. Computations were performed on the 128- 
node Pentium IV cluster at JLab, and at ORNL, 
under the auspices of the US DOE's SciDAC Ini- 
tiative. 
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Figure 1. The upper figure shows the effective 
masses for the three different embeddings of G\ g 
(smaller) and G\ u (larger) as the three different 
plotting symbols, and valence quark mass m = 
0.01. The points corresponding to the different 
embeddings are offset for clarity. The lower figure 
shows the effective masses for H u (diamond) and 
H g (square). 
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Figure 2. The upper figure shows the ra- 
tio (GigGi u ) I \G\ g G\ g ) , whilst the lower figures 

ShOWS (Gl^Gl™ 2 ) /{G\ g Wl G\ g Wl ). 



